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1 .  INTRODUCTION 

This  report  describes  the  initial  phase  of  a  program  which 
ultimately  will  allow  the  reliable  prediction  of  chemical  species  and 
temperature  within  the  atmosphere  following  a  low  altitude  nuclear 
explosion.  The  development  of  fireballs  from  low  yield  (<100  kton) 

bursts  at  these  altitudes  ( <30  km),  of  direct  concern  to  Site  Defense, 

will  be  dominated  by  buoyancy  forces.  The  initial  stage  of  these 
fireballs  resulting  from  large  energy  release  within  the  atmosphere  is 
characterized  by  spherical  symmetry,  the  domination  of  radiative  energy 
transport,  and  the  propagation  of  a  strong  spherical  shock  wave. 

Approximately  one  second  after  detonation,  the  fireball  comes  to  pressure 
equilibrium  with  the  ambient  air  and  because  of  its  low  density,  begins 
to  rise.  As  the  fireball  rises,  circulation  is  generated  and  a  toroid  is 
formed.  The  series  of  laboratory  experiments  currently  being  performed 

are  designed  to  simulate  these  "buoyant  rise"  fireballs  after  pressure 

equilibration. 

The  reliable  prediction  of  both  radar  and  optical  signatures 
resulting  from  atmospheric  nuclear  detonations  is  an  important  and  necessary 
prerequisite  for  a  defense  system  design  and  evaluation.  Before  such 
signatures  can  be  predicted,  the  resulting  spatial  distribution  of  temperature 
and  chemical  composition  must  be  known  and  this  requires  a  reliable 
theoretical  entrainment  model.  At  present,  no  data  exists  for  the  validation 
of  such  modeling.  Optical  data  from  nuclear  tests  give  only  gross  features 
of  the  events  such  as  fireball  size  and  rise,  but  what  is  required  is  a 
knowledge  of  the  temperature,  species,  and  velocity  fields.  This 
experimental  study  will  establish  a  sound  physical  data  base  to 
support  the  development  and  evaluation  of  theoretical  mixing  and  entrainment 
models. 
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2.  SUPPORT  ANALYSIS 

A  phenomenological  model  of  a  rising  fireball  or,  more  generally, 
an  isolated  buoyant  mass  has  been  developed.  The  main  emphasis  in 
development  of  this  model  was  oriented  toward,  although  not  limited  to, 
a  supporting  role  in  the  design  and  development  of  the  experimental 
fireball  simulation  facility.  For  this  reason  and  because  most  of  the 
physics  of  the  various  details  of  the  fireball  behavior  (such  as  the 
entrainment  of  external  fluid  and  subsequent  turbulent  mixing)  are 
not  well  understood,  the  approach  taken  was  to  satisfy  in  an  integrated 
sense  the  various  conservation  laws  of  mass,  momentum  and  energy. 

A  recent  paper  by  Wang  (Reference  1)  presented  an  analogous 
derivation  to  the  present  approach  as  have  Morton,  G.  I.  Taylor  and  Turner 
in  Reference  2  for  very  small  density  differences.  These  papers,  as 
important  and  informative  as  they  are,  are  not  necessarily  applicable  to 
either  the  fireball  or  the  laboratory  case  under  investigation  here. 

In  the  latter  analyses  the  effects  of  different  species  are  not  contained 
whereas  in  the  former  a  similar  situation  prevails  along  with  the  fact 
that  atmospheric  scale  height  effects  are  not  included. 

The  purpose  of  this  analysis  has  been  to  extend  these  works  to 
include  these  phenomena.  A  discussion  and  derivation  of  this  model  is 
given  below. 

2.1  GENERAL  REMARKS 

Consider  a  generally  spherical  region  with  density  formed  within 
an  external  atmosphere  with  a  larger  density  p^.  This  spherical  mass 
will  rise  initially  with  an  acceleration  of  2g's  because  of  the  coupled 
density  difference  and  gravitational  field.  Observations  (e.g..  Reference  3) 
indicate  that  this  acceleration  does  not  continue  indefinitely  because  of 
two  effects.  The  first  of  these  is  that  the  low  density  region  undergoes 
a  geometric  change  and  forms  into  a  torus  after  a  translation  of  approximately 
one  initial  diameter.  Commensurate  with  this  geometric  change,  the  buoyant 
region  begins  entraining  external  fluid  at  a  greatly  enhanced  rate  and 
thus  the  mass  of  fluid  carried  with  the  torus  increases  significantly  (the 
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buoyant  force  also  diminishes  for  displacements  on  tie  order  of  a  scale 
height).  It  is  this  phase  of  fireball  development  that  is  treated  here 
since  during  this  stage,  the  degree  of  turbulent  mixing  and  entrainment 
of  cool  air  into  the  fireball  will  establish  the  temperatures  and  chemical 
species  present  for  the  third  or  final  stage  of  the  fireball  history,  the 
"pancaking"  about  the  stabilization  altitude. 

After  the  initial  torus  formation  stage,  it  is  necessary  for 
mathematical  simplicity  to  assume  that  the  fireball  is  similarly  shaped 
during  the  period  of  interest.  The  fact  that  this  is  approximately  correct 
is  given  by  field  data  in  classified  reports  as  well  as  by  simple  experiments 
(see  Reference  4  for  instance).  From  these  latter  measurements,  it  is 
also  evident  that  the  fireball  (here  we  are  talking  about  the  entire 
region  carried  along  with  the  initial  products)  is  roughly  spherical  in 
shape.  Thus,  the  fireball  surface  area  and  volume  can  be  characterized 
by  a  single  parameter  'a'  which  can  be  thought  of  as  the  radius  of  a 
sphere  of  the  same  volume. 

2.2  EQUATIONS  OF  MOTION 

The  continuity  or  entrainment  equation  is  assumed  to  be  of  the  form 

dm. 

-fit  =  a  pm  u(4ira2)  (1) 

where  m  is  the  mass  of  entrained  atmospheric  air 

d 

pra(z)  is  the  atmospheric  density 

u  is  the  rise  velocity  of  the  fireball  center 

4ira2  is  the  surface  area  of  the  equivalent  sprtere 

a  is  the  dimensionless  entrainment  rate  which  is 
assumed  constant 

We  can  also  define  average  densities  pa  and  p.  from 

a  ! 
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/ 

m. 

and  Pi  =  (4/3)iraJ  2b 

where  m.  refers  to  the  initial  mass  of  the  fireball  and  is  a  given  constant. 
Obviously,  the  average  fireball  density  is  given  as 

p  =  p.  +  (3) 

a  I 

The  equation  of  motion  for  the  fireball  is 

|  ira3(p  +  \  p J  ^  =  j  ^(p*,  '  p)9  -  CD(^-  p^u2)  Ua2)  (4) 

The  term  ^  p^  within  the  parentheses  on  the  LHS  of  equation  (4)  accounts  for 
the  "apparent  mass"  effect  or  the  additional  inertia  caused  by  accelerating 
a  spherical  object  in  an  inviscid  fluid.  The  RHS  is  simply  the  buoyant  force 
less  the  drag  force.  The  drag  coefficient  CD  includes  both  the  form  drag 
CD'  and  the  entrainment  drag  given  by  8a.  It  is  assumed  that  Cp'  is 
constant.  This  assumption  is  not  critical  in  that  it  will  be  shown  that  the 
form  drag  is  small  compared  to  the  entrainment  drag. 

The  final  consideration  which  will  form  a  complete  system  is  the 
energy  equation.  To  formulate  this  to  the  same  degree  of  approximation  as 
the  previous  equations  it  is  assumed  that  the  pressure  within  the  fireball 
is  uniform  and  is  equal  to  the  ambient  pressure  at  the  center  elevation  of 
the  fireball.  This  assumption  does  become  weaker  for  fireball  dimensions 
on  the  order  of  an  atmospheric  scale  height  (^km) .  Further,  it  is  assumed 
that  the  kinetic  energy  is  small  compared  with  the  thermal  or  potential 
energy  of  the  fireball  (this  ratio  is  on  the  order  of  uVyRT^  which  is 
small).  The  total  enthalpy  within  the  fireball  is  then: 


Ht  =  |  »a3(p.Cp1  +  Pacpa)T 


where  c  .  and  c  are  the  specific  heats  of  the  initial  mass  and  the 
pi  pa 

entrained  air,  respectively. 
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For  moderate  pressures  and  temperatures,  the  perfect  gas  equation 
of  state  then  permits  this  to  be  written  as: 


H.  =  w-  Tra3p  -^r 
t  3  “  R 


(6) 


where  the  mass  concentrations  are  defined  as 


c  P<  m1 
p  m,  +  m 


•  c  -  —  -  — ^5— 

’a  p  m.  +  m 


(7) 


and  the  average  specific  heat  and  gas  constant  are,  respectively: 


c  =  C.c  .  +  C  c 
p  i  pi  a  pa 


R  =  C.R.  +  C  R 

1  a  a 


(8) 


The  rate  of  change  of  the  total  enthalpy  is  then  given  from  the  Second 
Law  of  Thermodynamics  and  assuming  that  the  entrainment  and  mixing  are 
adiabatic  over  the  entirety  of  the  fireball.  Thus, 


dH.  drn 

»  =  r  T  A  4. 

"eft"  -  cpa  a  dt 


(j  "a3> 


(9) 


The  second  term  on  the  RHS  accounts  for  the  fact  that  the  ambient  pressure 
varies  with  fireball  position  or  time.  This  variation  is  approximated  by 
assuming  that  locally  the  pressure-altitude  dependence  is  given  by: 


_L^P- .  1 
P.  dz  L 


(10) 


where  L  is  termed  the  pressure  scale  height  and  is  a  slowly  varying  function 
of  altitude  (e.g.,  L  ^  8.7km  at  sea  level  and  L  -v  6.9km  at  30km  altitude). 

Differentiating  equation  (6)  and  equating  it  to  equation  (9)  yields 
after  some  algebra 


da  1  au 
dt  "  "  3  yL 


auA 


(ID 
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where 


A  =  -P--  + 
cnRa  + 
P  a 


Ra/R  ~  cpa/cp 
p/p„ 


(12) 


Thus,  the  full  set  of  differential  equations  which  describe  the 
evolution  of  the  fireball  are: 


dm 

~df  =  “(4ira2)P<xu 


du  (p°° 

at =  — 


pa  -  Pj)9  -  f  cpP,u2/a 
(J  Pec  +  Pa  +  pi} 


da  1  au 

ar " ~  3 


+  auA 


(13) 


These  equations  can  be  easily  integrated  numerically  subject  to  the  initial 

conditions  at  t  =  t  that 

o 


mi 

3  =  ao*  u  =  V  pa  =  °»  pio  =  ( 4/3 ) rra^ 


(14) 


and  the  algebraic  relations 


pin(  a°)  and  Cn  =  cn'  +  8a 


(15) 


The  determination  of  the  initial  fireball  radius,  aQ  can  be 
ascertained  from  either  experimental  data  or,  in  the  case  of  a  nuclear 
explosion,  estimated  from  the  yield,  Y  (in  kilotons)  and  the  atmospheric 
density  at  burst  height,  Poo  (in  kg/m3).  An  example  of  such  an  expression 
taken  from  Reference  5  is: 
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The  "initial"  velocity  is  not  taken  as  zero  inasmuch  as  this 

.  0  i  ,  ;  i 

analysis  is  only  intended  for  times  after  torus  formation.  Although  the 
net  displacement  of  the  fireball  is  small  [0(ao)]  the  vertical  velocity  at 
this  location  is  not  negligible.  Investigations1 as  to  how  to  estimate  this 
are  still  underway.  Preliminary  indications,  however,  suggest  that  torus 
formation  occurs  when 


u„  ^  0.4ut 
o  i  T 


0.4 


I 


(16) 


where  uT  is  the  terminal  velocity  of  a  nonfentraining  sphere. 

i 

One  last  parameter  that  must  be  estimated  before  a  solution  can  bfe 
obtained  is  the  entrainment  constant  a.  Recourse  to  available  experimental 
data  is  necessary  to  determine  this  quantity  in  lieu  of  the  fact  that1  a 
theoretical  estimate  is  not  available  from  existing  turbulent  entrainment 
theories.  The  most  straightforward  estimate  of  a  can!  be  obtained  from  the 
experimental  data  for  the  fireball  size  as  a  function  of  position.  From 
the  last  of  equations  (13)  the  spatial  derivative  of  the  fireball  size  is 

I 

(17) 


If  only  one  species  is  involved,  the  quantity  A  is  identically  one 
(even  for  two  species  mixtures  it  is  found  that  the  entrained  mass  far 
overwhelms  the  initial  mass  so  that  this  approximation  is  also  valid).  Note 
that  in  equation  (17)  the  growth  rate  is  effected  by  scale  height  in  addition 
to  the  usual  entrainment  growth  rate. 


given  by 


da 

dz 


3yL 


+  aA 


where  the  fact  that  u 


dz 

at 


has  been  used. 


Data  from  small  scale  helium  release  experiments  in,  Reference  6  showed 
a  range  of  values  for  a  from  0.08  to  0.18  with  an  average  of  0.12.  Similar 
preliminary  experiments  were  performed  as  part  of  this  study  which  yielded 
values  for  a  from  0.10  to  0.15  with  an  average  value  of  0.125. 


18895-6002-R0-00 


Fielrf  test  data  from  nuclear  tests  show  typical  values  for  a  to 
vary  from  0.13  to  0.26  if  atmospheric  scale  height  effects  are  accounted 
for.  The  enormous  variation  in  scale  between  the  experiments  and  the  field 
data  should  not  be  overlooked  here--the  nuclear  explosions  are  some  100,000 
times  larger  than  the  experiments.  A  similar  variation  in  Reynolds  number 
i  also  noted  ( 1 0 **  for  the  lab  compared  to  10^  or  10^®  for  nuclear  explosions). 
In  spite  of  this,  the  small  variation  in  the  entrainment  constant  a  which 
is,  at  most,  a  factor  of  two  over  the  entire  range  of  Reynolds  number  is 
indeed  gratifying  although  the  size  effect  is  not  completely  removed. 

The  spherica  form  drag  coefficient  CD'  for  Reynolds  numbers  104  is 
roughly  on  the  order  of  0.2  or  lower.  Comparing  this  to  the  "entrainment 
drag"  coefficient  8a  shows  that  the  form  drag  is  less  than  10%  of  the  total 
drag  given  in  equation  (15).  Thus  the  variation  in  Cp'  with  Reynolds 
number  is  relatively  unimportant  insofar  as  the  dynamics  of  the  fireball 
is  concerned 

One  additional  feature  of  the  solution  to  these  equations  is  the 
prediction  of  the  stabilization  altitude.  This  is  the  altitude  where  the 
fireball  eventually  comes  to  rest  (actually  the  fireball  oscillates  about 
this  altitude  at  roughly  the  Brunt-Vasai la  frequency  oi  -  [-3-  -j^-]  ^  ).  Thus, 
the  solutions  are  valid  up  to  the  time  where  the  "collapse""of  the 
spherical  region  by  differential  buoyant  forces  becomes  significant.  For 
nuclear  events,  this  is  on  the  order  of  tens  of  minutes. 

The  equations  of  motion  of  the  fireball  derived  in  this  section 
have  been  programmed  for  numerical  integration  and  will  be  used  to  provide 
theoretical  results  for  comparison  with  both  laboratory  and  field  data. 
Additional  future  work  in  this  area  will  be  to  examine  these  data  to 
improve  estimates  of  the  initial  conditions  (t  ,  u  ,  a  ) 

O  O  0 

It  should  be  emphasized  again  that  the  primary  purpose  of  this  fireball 
analysis  is  to  support  the  experimental  study.  A  sophisticated  2-dinensional 
time  dependent  code,  well  beyond  the  scope  of  the  present  program,  would  be 
required  to  incorporate  the  turbulent  mixing  and  entrainm-nt  models  (as 
derived  from  the  experiments  for  example)  and  predict  the  distribution  of 
chemical  species  and  temperature  within  the  atmosphere  following  a  nuclear 
explosion. 
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3.  EXPERIMENTAL  FACILITY 


A  test  facility  consisting  of  several  primary  components  including 
a  high  pressure  test  tank  with  optical  ports,  a  pressurization  and 
exhaust  svstem  for  the  tank,  equipment  for  the  remote  formation  and 
controlled  bursting  of  helium-filled  soap  bubbles,  and  the  necessary 
diagnostic  equipment  to  obtain  density  and  velocity  field  measurements 
have  been  designed,  developed,  and  tested. 

3.1  HIGH  PRESSURE  TANK  AND  AUXILIARY  EQUIPMENT 

The  high  pressure  tank  is  an  ASME  pressure  vessel  with  a  diameter 
of  54  inches  and  a  length  of  14  feet.  A  photograph  of  the  test  tank 
installed  in  the  TRW  Fluid  Mechanics  Laboratory  is  shown  in  Figure  1. 

This  facility  was  designed  to  operate  in  a  pressure  range  from  vacuum 
conditions  to  10  atmospheres. 

Viewing  ports  were  designed  to  provide  a  12  inch  diameter 
horizontal  view  and  a  10  inch  diameter  vertical  view.  The  horizontal 
view  ports  are  located  10  inches  above  the  tank  centerline  to  allow  torus 
vertical  trajectories  of  15  bubble  diameters  (based  on  an  initial  soap 
bubble  diameter  of  2  inches)  without  significant  wall  interference. 

The  high  pressure  tank  is  equipped  with  both  a  pressurization  and 
venting  system.  Tank  venting  is  accomplished  through  a  symmetrical  duct 
system  consisting  of  4  inch  diameter  flanges  on  each  end  of  the  tank 
which  are  piped  to  a  common  high-speed  valve  and  muffler  system  outside 
the  building.  The  symmetrical  vent  system  and  control  valve  are  shown  in 
Figure  1.  The  high  pressure  tank  was  designed  symmetrically  and  with  a 
length-to-diameter  ratio  of  about  3  to  minimize  induced  velocities  on 
the  torus  motion  during  tank  venting.  The  vent  crntrol  valve  is  a  6  inch 
diameter,  pneumatically  actuated,  Lunkenheimer  ball  valve  which  will 
allow  tne  tank  pressure  to  be  reduced  from  10  to  2  atmospheres  in  about 
2  seconds.  High  pressure  nitrogen  bottles  manifolded  together  are  used 
to  pressurize  the  tank. 
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A  system  was  designed  and  developed  to  remotely  form  and  burst 
helium-filled  soap  bubbles  in  the  high  pressure  tank.  Positive  displacement 
pistons  and  a  valving  system  are  used  to  meter  the  amount  of  helium  and 
soap  solution  to  the  bubble  tube  to  accurately  control  the  bubble  size. 

A  small  electrode  spaced  about  20  mils  from  the  bubble  tube  has  been  used 
with  a  200  microfarad  capacitor  discharge  system  to  repeatedly  burst  the 
bubbles.  Several  other  methods  of  bursting  the  bubbles  are  currently 
being  investigated  to  Improve  if  possible  the  symmetry  of  the  Initial 
helium  release.  The  bubble  formation  head  and  present:  electrode  breaker 
have  been  designed  as  an  Integral  unit  so  they  can  be  traversed  vertically 
within  the  tank  through  a  1  inch  diameter  access  flange  and  seal.  The 
bubble  tube  position  can  be  varied  from  the  center  to  30  inches  below  the 
horizontal  viewing  ports. 

3.2  DIAGNOSTIC  EQUIPMENT 

In  order  to  completely  specify  the  fluid  state  and  motion  the 
following  must  be  known  for  a  non- reacting  gas:  two  thermodynamic 
variables,  such  as  pressure  and  temperature;  a  variable  of  the  motion,  such 
as  velocity;  and  the  composition  of  the  field.  Measurements  of  the  tank 
pressure  and  temperature,  outside  the  induced  flow,  suffice  for  specifying 
the  pressure  and  temperature  of  the  field  since  the  Induced  fluid  motion 
Is  Incompressible  (1.®.,  low  velocity  and  adlabaMc).  The  remaining 
unknowns,  density  (or  composition)  and  velocity  are  measured  optically 
using  holographic  Interferometry  and  a  particle  tracking  technique, 
respectively.  Optical  techniques  were  adopted  for  this  study  because  of 
three  primary  difficulties  Involved  with  flow  field  probes.  First, 
Interpretation  of  data  obtained  from  probes  In  a  highly  three-dimensional 
flow  Is  almost  Impossible.  Also,  probe  interference  effects  are  sizable 
in  an  Incompressible  flow  such  as  this,  and  finally  the  cost  and  time 
consumed  In  reconstructing  a  time  dependent  flow  field  from  local 
measurements  are  prohibitive. 
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3.2.1  Holographic  Interferometer 

Pulsed  laser  holographic  Interferometry  will  be  used  to  measure 
the  mean  density  field  and  turbulent  structure  of  the  rising  toroids. 

This  technique  has  been  used  extensively  at  TRW  to  measure  mean  and 
turbulent  flow  field  properties  of  wakes  as  described  In  Reference  7.  The 
laser  holograph*  which  will  be  used  In  this  program  Is  similar  In  design 
to  the  one  used  In  previous  studies  but  has  been  considerably  Increased 
In  size  to  accommodate  the  large  high  pressure  tank.  The  main  feature  of 
the  holocamera  Is  that  It  records  orthogonal  views  simultaneously  across 
the  tank  test  section.  Consequently,  this  system  Is  capable  of  recording 
any  transmission  subject  In  a  9  Inch  diameter  area  In  line  with  the  8  by 
10  Inch  photographic  plates  and  the  large  18  Inch  objective  lens  at  the 
opposite  ports.  The  advantages  of  using  this  holographic  interferometer 
as  compared  to  the  standard  Mach  Zehnder  Interferometer  include  the  three- 
dimensional  viewing  and  depth  focusing  properties  of  holographic  Interfero 
grams,  and  the  fact  that  precision  optics  and  accurate  alignment  are  not 
required  with  holograph;:  Interferometry  since  a  common  optical  path  Is 
used  for  both  the  scene  and  comparison  beams.  An  end  view  of  the  tank  and 
ports  and  holocamera  components  In  close  proximity  to  the  tank  are  shown 
In  Figure  2  while  a  schematic  of  the  complete  holocamera  Identifying  all 
primary  optical  components  Is  presented  In  Figure  3.  A  photograph  of  the 
holocamera  together  \  1th  the  test  tank  Is  shown  In  Figure  4. 

In  this  holocamera  the  orlmary  beam  for  the  laser  Is  split  into 
a  reference  and  a  scene  beam.  The  scene  beam  Is  directed  through  the 
scene  to  the  plates  whereas  the  reference  beam  Is  directed  on  the  plates 
without  passing  through  the  disturbance.  The  prism  plate  and  lenses  shown 
In  Figure  3  are  used  to  spatially  and  temporally  match  the  scene  and 
reference  beams  over  the  plates.  A  Q-swItched  solid  state  ruby  laser  with 
an  energy  of  about  1  joule  and  pulse  length  of  about  100  nanoseconds  Is 
used  to  expose  both  the  vertical  and  horizontal  beam  holograms  which  are 
recorded  on  Agfa-Gavaert  8E75  plates. 


*The  term  laser  holograph  pertains  to  both  the  laser  cavity  or  oscillator, 
and  all  other  optical  components  which  comprise  the  holocamera. 
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Interferograms  are  obtained  in  this  common  path  interferometer  by 
sequentially  recording  the  comparison  and  scene  beams  on  the  same  plates. 
The  double  exposure  of  the  plates  by  the  two  pulses,  the  one  of  the 
undisturbed  volume  and  the  other  of  the  same  volume  with  the  test  scene 
introduced,  produces  the  interferograms.  The  standard  mode  of  operation  of 
this  holocamera  in  which  the  position  of  all  optical  components  remain 
fixed  during  the  time  between  the  two  exposures  produces  an  Interference 
record  called  an  Infinite  fringe  Interferogram.  A  very  small,  accurately 
controlled  rotation  of  a  mirror  and  the  corresponding  plate  between  the 
two  exposures,  will  yield  a  different  type  of  Interference  pattern  on 
the  plate  called  a  finite  fringe  Interferogram.  The  hardware  necessary 
to  permit  this  accurate  rotation  of  optical  components  has  been  designed 
and  fabricated  for  the  holocamera  horlzonta.  beam  and  will  be  used 
throughout  the  study  as  necessary.  A  discussion  of  the  interpretation  of 
both  finite  and  infinite  fringe  Interferograms  and  the  advantages  and 
disadvantages  of  each  Is  included  In  Section  4.2. 

3.2.2  Particle  Tracking  Apparatus 

The  fluid  velocity  will  be  measured  by  photographically  recording 
the  motion  of  tracer  particles  which  follow  the  fluid  motion.  Provided 
that  the  particle  Reynolds  nimber  based  on  Tag  velocity  and  particle 
diameter  Is  less  than  about  1,  the  drag  force  is  viscous  rather  than 
inertia  dominated,  and  the  particles  will  quickly  attain  and  follcw  the 
fluid  motion.  An  exhaustive  treatment  of  using  particulate  tracers  for 
fluid  velocity  measurements  was  conducted  In  a  wind  tunnel  by  Bourot  as 
described  In  Reference  8.  More  recently,  flow  visualization  studies  have 
successfully  been  conducted  utilizing  light  reflection  from  small  helium- 
filled  soap  bubbles  as  Indicated  In  Reference  9.  During  the  course  of 
this  recent  work,  equipment  was  developed  by  Sage  Action,  Inc.  (SAI)  which 
could  continuously  produce  helium-filled  soap  bubbles  as  tracer  particles 
with  various  degrees  of  buoyancy. 

An  SAI  Model  3  console  and  low  speed  bubble  head  will  be  used  In 
the  current  program  to  produce  slightly  buoyant  bubbles  about  1/8  Inch  In 
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diameter.  The  bubble  head  will  be  mounted  at  the  bottom  of  the  test  tank 
on  a  traversing  mechanism  so  a  vertical  sheet  of  bubbles  about  lft  inches 
long  will  be  produced  along  the  tank  axis. 

A  1500  watt  tungsten-halogen  lamp  rated  at  30,000  lumens  will  be 
attached  above  the  top  view  port  and  masked  to  produce  a  high  intensity 
light  beam  1/4  Inch  wide  by  about  10  Inches  long  coincident  with  the  sheet 
of  bubbles.  A  4  by  5  Inch  view  camera  will  be  mounted  at  one  of  the 
side  windows  and  focused  on  the  light  sheet  to  record  the  motion  of  all 
Illuminated  bubbles  In  the  field  of  view.  A  *  4 "1 mum  of  three  multiple 
exposures  will  be  obtained  In  each  test  run  with  an  exposure  time  of 
about  100  msec  using  a  MAT I  electroshutter.  Each  film  record  will  show 
a  particle  trajectory  as  a  row  of  dashes  the  lengths  of  which  will  be 
used  to  determine  the  local  velocity  at  that  point  In  the  field.  Velocity 
data  will  only  be  reduced  from  central  dashes  In  each  particle  trajectory 
since  the  Initial  and  final  dash  In  a  trajectory  could  be  produced  with 
a  particle  either  entering  or  leaving  the  light  beam  during  the  exposure. 
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4.  DATA  ACQUISITION  AND  INTERPRETATION 


The  purpose  of  this  section  Is  to  discuss  both  the  procedure  for 
obtaining  holoaraphlc  Interferograms  and  the  physical  Interpretation  of 
this  form  of  data.  Acquisition  of  velocity  field  data  with  the  particle 
tracking  apparatus  was  discussed  previously  In  Section  3.2.2. 


4.1  ACQUISITION  OF  INTERFEROMETRY  DATA 

The  procedure  which  Is  used  to  obtain  an  orthogonal  set  of  Infinite 
fringe  Interferograms  of  a  simulated  fireball  consists  of  the  following 
primary  steps: 

1.  Plate  holders  containing  the  special  fine  grain,  low 
speed  photographic  plates  arc  loaded  In  the  horizontal 
and  vertical  plate  holder  receptacles. 

2.  The  slides  covering  the  plates  are  drawn  after  the  laser 
and  bubble  blowing  apparatus  have  been  checked  and  all 
sources  of  Illumination  In  the  test  area  have  been 
eliminated. 

3.  A  helium-filled  soap  bubble  Is  blown  to  a  predetermined 
size  and  the  laser  power  supply  Is  charged. 

4.  About  five  seconds  after  Its  formation,  the  helium 
bubble  Is  burst. 

5.  Simultaneous  with  the  bubble  bursting,  an  electrical 
signal,  first  sent  to  a  time  delay  circuit,  activates 
the  laser  firing  mechanism  which  records  the  test 
scene  on  the  plates.  The  time  delay  Is  adjusted  to 
allow  the  helium  cloud  to  rise  to  near  the  center  of 
the  windows.  A  typical  time  delay  of  650  milliseconds 
would  allow  a  bubble  with  an  Initial  diameter  of  one 
Inch  to  rise  about  five  Inches  at  a  tank  pressure  of 
eight  atmospheres. 

6.  A  comparison  scene  Is  recorded  on  the  plates  by  firing 
the  laser  a  second  time  about  30  seconds  after  the 
event  scene.  This  delay  ensures  that  the  disturbance 
Is  well  outside  of  the  viewing  region  resulting  In  a 
uniform,  constant  density  comparison  scene. 
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7.  The  photographic  plates  are  removed  and  developed  using 
special  photographic  developing  agents. 

8.  The  holograms  are  reconstructed,  l.e.,  observed  and 
photographed  by  Illuminating  the  holograms  with  a  cw 
laser. 

9.  Large  scale  working  photographs  are  now  produced  for 
use  In  the  data  reduction  phase  as  described  In 
Section  5.1.3. 

A  horizontal  view  finite  fringe  Interferogram  Is  obtained  by 
following  the  exact  procedure  as  outlined  above  except  for  the  addition  of 
one  step.  This  consists  of  rotating  the  horizontal  beam  plate  holder  and 
adjacent  mirror  a  small  precise  amount  between  steps  5  and  6  above,  l.e., 
between  recording  the  test  and  the  comparison  scenes  on  the  plate. 

4.2  INTERPRETATION  OF  HOLOGRAPHIC  INTERFEROGRAMS 

Interpretation  of  both  Infinite  fringe  and  finite  fringe  Intcrferograms 
will  be  Included  In  this  section  since  both  types  of  data  will  be  acquired 
and  analyzed  during  the  course  of  this  study. 

The  principle  of  operation  of  an  Interferometer  Is  that  optical 
Interference  occurs  between  two  coherent  beams  passing  through  media  of 
different  density.  Whenever  optical  paths  are  the  same  or  differ  by  an 
Integral  multiple  of  X,  (where  x  Is  the  wavelength  of  the  monochromatic 
light  source.  6943A  In  this  study)  the  beams  reinforce  one  another  whereas 
when  a  disturbance  which  shifts  the  optical  path  by  any  multiple  of  x/2  Is 
encountered,  complete  Interference  occurs.  Complete  reinforcement  will 
show  up  as  light  fringes  on  a  hologram  plate  and  complete  Interference 
as  dark  fringes. 

4.2.1  Infinite  Fringe  Interferoqrams 

An  example  of  an  Infinite  fringe  Interferogram  Is  presented  In 
Figure  5.  This  horizontal  view  hologram  was  obtained  with  an  unbroken 
one  inch  diameter  helium  bubble  at  a  tank  pressure  of  one  atmosphere.  Each 
fringe  In  this  type  of  Interferogram  corresponds  to  a  locus  of  constant 
optical  path  so  Interpretation  for  data  reduction  purposes  consists  of 
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numbering  and  locating  the  fringes  along  any  plane  of  interest  through 
the  fireball.  These  fringe  numbers  and  locations  can  directly  be  used  to 
determine  the  integrated  density  or  helium  concentration  change  through 
the  fireball  along  that  optical  path.  Details  of  the  computational 
procedure  required  to  obtain  concentration  profiles  from  these  fringe 
number  data  are  presented  in  Section  5.1.3.  In  the  sample  infinite  fringe 
interferogram  presented,  the  density  is  constant  inside  the  bubble  and 
lower  than  the  constant  density  region  outside  the  bubble  by  about  a 
factor  of  7.  Consequently,  the  interference  fringes  which  consist  of 
concentric  circles  with  decreasing  spacing  as  their  radii  increase  depict 
only  changes  in  optical  path  length  through  the  helium  sphere. 

A  second  example  of  an  infinite  fringe  interferogram  is  shown  in 
Figure  7.  This  interferogram  shows  a  vortex  formed  from  an  initial  one 
inch  diameter  helium  bubble  at  one  atmosphere  pressure  after  a  rise  of 
about  4  initial  bubble  diameters.  The  fringe  pattern  in  this  figure  is 
more  typical  of  that  expected  in  this  study  in  that  the  fringe  geometry 
and  spacing  result  from  a  toroidal  cloud  containing  a  mixture  of  helium 
and  air. 

One  disadvantage  of  Infinite  fringe  interferograms  is  that 
ambiguities  can  arise  In  reading  the  fringe  numbers,  particularly  if  the 
test  event  Is  turbulent.  Since  infinite  fringe  Interferograms  indicate 
the  number  of  wavelengths  that  emitted  light  advances  or  is  retarded  in 
phase  as  it  passes  through  the  disturbance,  it  may  not  be  apparent  whether 
a  fringe  indicates  an  Increase,  decrease,  or  no  change  in  fringe  number. 
Unless  there  is  an  a  priori  knowledge  of  the  flow  field  (obviously  there 
is  very  little  when  the  fireball  is  turbulent),  ambiguities  may  occur. 

This  fringe  ambiguity  problem  can  be  eliminated  by  recording  finite 
fringe  Interferograms  which  are  described  in  the  following  section. 

4.2.2  Finite  Fringe  Interferograms 

An  example  of  a  finite  fringe  interferogram  is  presented  in 
Figure  8.  This  Interferogram  was  obtained  at  the  same  test  conditions  as 
the  infinite  fringe  interferogram  in  Figure  7.  In  this  type  of  hologram, 
the  test  disturbance  causes  a  displacement  or  shift  In  the  fringes  relative 
to  the  undisturbed  portion  of  the  test  section.  An  increase  in  optical 
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path  Is  shifted  In  the  opposite  direction  to  a  decrease  In  optical  path, 
consequently  no  ambiguity  ca*'  occur  In  regard  to  the  fringe  shift  direction. 
The  magnitude  of  the  fringe  shift  In  this  type  of  hologram  can  be 
directly  related  to  fringe  number  in  an  infinite  fringe  hologram  and 
therefore  to  an  Integrated  change  In  density  or  concentration  along  the 
optical  path. 

A  second  example  of  a  finite  fringe  interferogram  is  presented  in 
Figure  6.  This  hologram  of  another  unbroken  helium  bubble  was  obtained 
at  the  same  conditions  as  the  infinite  fringe  interferogram  of  Figure  5. 

The  spacing  of  the  parallel  fringes  in  this  hologram  is  too  large  to 
allow  the  fringes  to  be  just  shifted  as  they  were  in  Figure  7.  Instead, 
several  elliptic  shaped  fringes  were  formed  near  the  top  of  the  bubble. 
Although  it  Is  more  difficult,  even  these  fringes  can  be  interpreted  and 
the  data  reduced  since  the  shift  of  each  fringe  will  be  relative  to  a 
specific  determinable  undisturbed  fringe  above  the  bubble.  Finite  fringe 
holograms  with  easier  to  reduce  continuous  fringes  can  be  obtained  in 
small  events  with  high  helium  concentrations  (such  as  this  unbroken  bubble) 
by  increasing  the  amount  of  mirror  and  plate  rotation  between  recording  the 
test  scene  and  comparison  scene  on  the  hologram  plates. 

One  disadvantage  of  this  form  of  hologram  is  that  a  visual  picture 
including  the  outline  of  the  test  event  is  not  obtained,  however  this  is 
more  than  offset  by  the  simplifications  in  the  data  reduction. 
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5.  DATA  REDUCTION 

The  purpose  of  this  section  is  to  describe  briefly  the  technique 
for  reducing  holographic  interferometry  data  as  it  applies  to  the  present 
study.  The  process  of  achieving  the  final  result,  namely,  determining  the 
mean  and  fluctuating  density  field  from  the  recorded  holographic 
interferograms  is  comprised  of  two  primary  steps:  first  fringe  interpretation, 
numbering  and  reading  from  the  plates,  and  then  computer  data  reduction 
of  the  density  field.  This  approach  has  been  used  in  a  number  of  former 
studies  (References  10,  11,  and  12)  which  involved  data  reduction  of 
both  mean  and  fluctuating  (turbulent)  model  reentry  vehicle  flow  fields. 

In  each  of  the  previous  studies,  density  changes  occurred  as  a 
result  of  pressure  and  temperature  changes.  In  the  present  study, 
density  changes  result  because  of  changes  in  molecular  weight  which 
occur  when  diffusion  and  mixing  take  place  between  helium  and  air  at 
constant  temperature  and  pressure.  The  interferometry  measurements  are 
basically  measurements  of  the  changes  in  the  index  of  refraction  field  which 
is  expressed  in  terms  of  both  the  mass  fraction  of  helium,  C^,  or  air, 

C^,  and  mass  density,  p.  The  perfect  gas  law  can  be  expressed  in  terms 
of  the  same  two  unknowns  and  thus  a  solution  for  the  species  mass 
fraction  and  density  can  be  obtained.  The  details  to  this  approach  for 
mean  quantities  is  given  in  Section  5.1  while  a  discussion  of  the  procedure 
for  obtaining  turbulent  data  is  included  in  Section  5.2. 

5.1  CALCULATION  OF  DENSITY  PROFILES 

5.1.1  Interfer.  letry  Theory 

The  standard  equation  for  fringe  shift  is 
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where  S  =  fringe  number 

X  =  wave  length  of  light  in  vacuum 
L  =  integration  path  length 
n  =  index  of  refraction 
s  =  path  length  along  the  light  ray 

subscript 

~  =  undisturbed  conditions  recorded  in  comparison  scene  and  also 
observed  in  „he  test  scene  as  a  reference  condition 

The  Abel  inversion  integral  (discussed  later  in  Section  5.1.2)  is 

used  to  solve  for  n  as  a  function  of  S  for  the  case  of  axial  symmetry. 

Here,  one  wants  to  consider  how  the  solution  n  -  n  can  be  used  to 

00 

calculate  concentration  or  mass  fraction  profiles  (r)  from  the 
interferograms. 

For  a  multi -component  gas  case,  the  index  of  refraction  becomes, 

n  =  E  V,  +  '  = »  E  ♦  >  .  pK  ♦  i 

where  K  =  Gladstone-Dale  constant 
p  =  density 
subscript 

i  =  species  i 
superscript 

=  mean  quantity 

For  helium  and  air  (subscripts  H  and  A,  respectively) 

n  -  1  =  P  (KaCa  +  KhCh> 

and  because 

Ec,  ■  i 


p  K 

rnn  r 


PKA  kh 

A  [CH  f-  -  1)  +  1]  -1 


P  K 


(1) 
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From  the  perfect  gas  law  the  density  ratio,  p/p^,  may  be  obtained  as 


P  p  T  M 

co  r  oo  'oo 


where 

T  =  temperature 
p  =  pressu.e 
M  =  molecular  weight 

A  simplifying  assumption,  which  in  any  event  is  non-crucial  to  the  final 
results,  is  that  p  =  p^  and  T  =  T^,  i.e.,  the  pressure  and  temperature  of 
the  comparison  and  test  scenes  are  the  same.  Then,  p/Poo  =  M/M^. 

From  Dalton's  Law  of  partial  pressures, 

C. 

P  «  TDPi  R1  =  pT/?£  Mi=  pT^~ 


where 

R  =  gas  constant  per  unit  mass 
R  =  gas  constant  per  mole 
and  M  is  defined  as 


For  hel ium  and  air, 


Consequently  the  density  ratio  becomes 
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The  case  here  Is  for  a  comparison  scene  (sub  °°)  comprised  of  air,  thus 

I 

M  =  M.  and  K  =  K. 
oo  a  °°  A 


and 


P_  _ 
p„  " 


V 


'H 


(H 


(2) 


+  1 


Eliminating  p/p^  from  equations  (1)  and  (2)  on®  obtains,  after  some 
algebra,  *  . 


-  n) 
pI  K 


CH  =  r(n  -  n  ) 


fTn-HooT  1  /"A  \  /?h  A 

RPC-  J  '  / '  W  ■  V 


(3) 


The  mole  fraction  of  helium  becomes: 


i  C  - 

lh  Ml 


xh  =  -K 


(4) 


where 


M 


i  i 


+  1 


and  the  density  ratio  becomes 


p  M 

=  ma 


(5) 
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The  quantity  Kji^  In  Equation  3  is  linearly  dependent  on  pressure, 
and  only  slightly  dependent  on  temperature  and  the  wavelength  of  light 
used.  A  formula  established  by  Biot  and  Arago  expressing  these 
dependencies  and  confirmed  by  more  recent  experiments  (Reference  13)  is 


where  n„  is  the  index  of  refraction  at  T  =  0°C,  p  =  1  atm 

0  oo  'oo 

a  is  the  coefficient  of  thermal  expansion  in  C0”1 
is  the  ambient  temperature  in  Cc 

Pot  is  the  ambient  pressure  in  atmospheres 


Values  of  nQ  and  a  required  for  this  study  were  obtained  from  tables  in 
Reference  13. 


5.1.2  Abel  Inversion  and  Numerical  Integration 


The  basic  interferometry  equation  which  relates  fringe  number  or 
fringe  shift  to  the  integrated  change  in  index  of  refraction  through  a  flow 
field  is 


where  all  auantities  were  defined  in  the  previous  section.  Restricting 
consideration  to  an  axisymmetric  flow,  the  equation  becomes 


The  new  variables  are  defined  in  the  figure  below: 
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To  solve  for  n  -  n  ,  one  notes  that  this  equation  Is  the  same  as 

oo 

Volterra's  Integral  equation  of  the  first  kind,  with  the  exception  that 
the  limits  of  Integration  are  different.  The  details  of  the  solution 
(Reference  14)  are  rather  lengthy  and  will  be  omitted  here.  It  becomes 


I"  -  Or-y 
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y  (r2  -  y2) 

One  can  express  this  Integral  as  a  sum,  ^Reference  15), 


(n  -  n  ) 


(S*  ~  S«tl)  {rK  +  M2  . 
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.  21  (W 

\lT 
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d(r2  -  y2) 1/2 


U2  -  12)1/2  (6) 


where  y^  =  1a 

t  =  zone  width 

i  =  0,  1,  2  ...  N  -  1  (zones) 

The  calculation  geometry  Is  centered  on  the  axis  of  symmetry  and 
Is  comprised  of  N  zones  having  constant  properties  as  shown  In  the  sketch: 


r  radius  to 
outerzone 


y,  boundary  of  1 
zone  where  n  Is 
evaluated 


l  calculation 
zone  size 


th 
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Having  n  -  n^  It  Is  a  simple  matter  with  a  computer  program  to  calculate 
helium  mole  fraction,  helium  concentration,  and  density  distributions 
using  Equations  3,  4,  and  5,  respectively  from  Section  5.1.1. 

5.1.3  Computational  Procedure 

Several  steps  are  required  In  the  data  reduction  procedure  to 
obtain  fringe  number  or  fringe  shift  distributions  necessary  for  Input  to 
the  Abel  Inversion  computer  program. 

Photographs  are  obtained  of  the  original  8  by  10  Inch  holographic 
plates  (both  horizontal  and  vertical  views)  yielding  4  by  5  Inch  negatives. 
These  negatives  are  enlarged  by  a  factor  of  approximately  ten  resulting 
In  large  scale  positive  transparencies.  A  horizontal  plane  coincident 
with  the  maximum  diameter  of  each  vortex  Is  selected  and  Identified  by 
scribing  a  line  on  the  working  photographs.  The  light  and  dark  fringes 
are  located  along  each  line  and  recorded.  Fringes  are  then  numbered 
starting  with  1/2  for  the  outermost  dark  fringe  (corresponding  to  the 
half  wavelength  shift  In  optical  path  discussed  In  Section  4.2),  1 
for  the  next  Inward  light  fringe  and  so  on,  Increasing  by  1/2  each  time 
a  new  fringe  is  encountered.*  Next,  the  location  of  the  vortex  center  Is 
determined  by  the  position  of  the  minimum  Internal  fringe  number.  These 
fringe  nunber  distributions  for  each  half-vortex  (having  Its  origin  at  the 
vortex  center)  provide  the  Input  to  a  computer  program  which  computes 
density,  helium  mole  fraction,  and  helium  concentration  as  functions  of 
radial  position. 

To  date,  all  data  reduction  has  been  performed  following  the 
procedure  outlined  above.  The  feasibility  of  automating  the  data 
reduction  by  digitizing  each  Inter ferogram  and  analyzing  the  digital  data 
with  a  sophisticated  computer  program  Is  being  studied.  The  DNA  Film 
Digitizing  Equipment  at  AFWL  operated  by  ISI  (see  Reference  16)  could 

♦Since  many  of  the  fringes  appear  as  multiply  connected  regions,  the  same 
fringe  may  Intersect  the  plane  In  question  more  than  once.  When  this 
occurs,  care  must  be  taken  to  assign  the  same  number  to  the  same  fringe. 
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possibly  be  used  for  this  purpose.  The  primary  difficulties  at  the 
present  time  are  the  software  developments  necessary,  both  In 
digitizing  and  analyzing  the  data. 

5.2  CALCULATION  OF  1URBULENT  STATISTICS 

Oata  reduction  of  turbulent  flows  from  holographic  Interferograms 
has  been  demonstrated  for  wakes  In  ballistic  range  firings  (References  10 
and  11).  Various  turbulent  statistical  quantities  were  calculated  using 
a  data  reduction  procedure  developed  especially  for  Interferograms  of 
hypersonic  wakes  (Reference  17).  The  main  feacure  of  this  technique  Is 
that  a  cross  view  localization  of  the  turbulent  statistics  Is  possible. 

Either  by  obtaining  data  at  two  angles  from  one  hologram  or  by  simultaneously 
obtaining  data  from  two  holograms  recorded  at  different  angles,  these 
data,  when  cross  correlated,  can  be  shown  to  provide  turbulent  statistics 
localized  In  the  vicinity  of  the  Intersection  of  the  two  views.  It  Is 
anticipated  that  with  a  straightforward  extension  of  this  technique  that 
at  least  some  of  the  turbulent  statistics  of  the  rising  fireball  can  be 
determined. 

In  this  procedure,  the  holograms  are  first  photographed  from  the 
desired  view  angle,  the  reconstructions  accomplished  using  a  cw  laser. 

From  the  negatives,  patterns  of  fluctuating  optical  density  (Indicative 
of  change  In  fringe)  are  recorded  by  a  microdensitometer.  The  extraction 
of  the  mean  flow  provides  a  new  aspect  to  the  data  reduction  due  to  the 
different  structure  of  the  fireball  phenomena  as  compared  to  the  turbulent 
wake.  Once  this  operation  has  been  accomplished,  only  fluctuations  due 
to  turbulence  remain.  These  turbulent  fluctuations  provide  the  Input  to 
a  fast  Fourier  transform  computer  program  which  was  developed  during  the 
aforementioned  wake  studies.  This  program  compiles  both  the  autospectrum 
and  autocorrelation  for  each  view,  and  the  cross  spectrum  and  cross 
correlation  as  well  as  an  Integral  length  scale  for  each  pair  of  cross 
views. 


In  the  future  an  attempt  will  be  made  to  also  obtain  turbulent 
velocity  fluctuation  data  using  the  particle  tracking  apparatus.  Again  the 
same  approach  must  be  used,  namely  extracting  the  mean  velocity  field  from 
the  data  to  yield  only  the  turbulent  components. 
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6.  EXPERIMENTAL  RESULTS 


A  discussion  of  both  qualitative  results  from  hardware  development 
tests  and  reduced  data  from  holographic  Inter ferograms  obtained  at  both  4 
and  8  atmospheres  are  presented  In  the  sections  below. 


6.1  FACILITY  DEVELOPMENT  AND  CHECKOUT  EXPERIMENTS 

Several  preliminary  tests  were  performed  to  support  the  development 
of  the  test  facility,  auxiliary  hardware  and  diagnostic  equipment.  Fireball 
"size  and  rise"  data  obtained  from  shadowgraph  motion  pictures  were  used 
to  establlsn  the  size  of  the  test  tank  and  optical  ports.  These  tests  were 
performed  at  one  atmosphere  resulting  In  entrainment  coefficients,  (see 
Section  2.2),  from  0.10  to  0.15  with  an  average  value  of  0.125*. 

Additional  shadowgraph  data  were  obtained  In  the  controlled  test  tank 
environment  during  the  checkout  phase  of  this  equipment.  A  laser  shadowgraph 
of  a  laminar  vortex  formed  by  the  release  of  helium  from  a  2  Inch  diameter 
bubble  at  one  atmosphere  Is  presented  In  Figure  9.  The  Reynolds  nurr^er** 
for  this  event  Is  3700.  It  Is  apparent  from  Figure  9  that  the  physics  essential 
to  toroid  formation  are  contained  within  this  experimental  simulation.  A 
turbulent  vortex  obtained  with  an  Initial  helium  bubble  diameter  of  1.1 
Inches  and  a  tank  pressure  of  8  atmospheres  corresponding  to  a  Reynolds 
number  of  12,300  Is  shown  In  Figure  10.  The  higher  Reynolds  number  vortex 
In  Figure  10  clearly  shows  both  a  large  and  small  scale  turbulent  structure 
whereas  the  vortex  In  Figure  9  obtained  at  a  lower  Reynolds  number  Is 
evidently  laminar  (l.e.,  appears  to  totally  lack  such  structure).  The 
evolution  of  a  vortex  from  a  sphere  of  helium  at  a  high  Reynolds  number  Is 
shown  as  sequential  frames  of  a  shadowgraph  motion  picture  film  In  Figure  11. 

It  Is  Interesting  to  note  that  structure  Identifiable  on  the  periphery  of 


*  These  preliminary  tests  were  performed  In  an  uncontrolled  environment  (l.e., 
subject  to  air  currents  In  a  laboratory)  which  may  account  for  some  of 
the  variation  In  the  fireball  growth  rate. 


♦♦Reynolds  nunber  Is  defined  as  Re  ■  p.u 
velocity  of  a  non-entraining  sphere,  0 
and  pm  and  u„  are  the  ambient  air  dens 
respectively. 


Dj/w.  where  uj  Is  the  terminal 
is  the  Initial  bubble  diameter, 
ty  and  dynamic  viscosity. 
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the  event  is  swept  from  the  top  around  to  the  bottom  and  is  either 
carried  Into  the  vortex  because  of  circulation  or  trails  to  form  a 
wake. 


In  addition  to  the  preliminary  experiments  just  described,  numerous 
bench  tests  were  conducted  to  Investigate  various  remote  soap  bubble 
forma tlor  and  bursting  techniques  which  resulted  In  the  development  of 
the  equipment  described  earlier  In  Section  3.1.  Recent  tests  utilizing 
Interferometry  data  obtained  Immediately  after  bubble  breaking  Indicate 
that  the  capacitor  discharge  bursting  system  Introduced  Initial  asyrmetrles 
In  the  hellun  release  which  were  not  detected  In  the  more  qualitative 
shadowgraph  data.  A  discussion  of  the  possible  effects  of  these  asynnetrles 
on  the  preliminary  data  Is  Included  In  Section  6.2  while  a  recently 
developed  Initialization  system  which  yields  a  symmetric  helium  release  Is 
described  In  Section  7. 

6.2  HOLOGRAPHIC  INTERFEROMEO  DATA 

It  should  be  noted  that  all  preliminary  fringe  distribution  and 
helium  concentration  results  presented  In  this  section  are  based  on 
Interferometry  data  obtained  In  tests  In  which  the  helium  release  was 
Initiated  with  the  electrode  system  described  in  Section  3.1.  It  Is 
expected  that  the  asymmetries  attributed  to  this  type  of  hellixn  release  will 
be  considerably  reduced  through  the  use  of  the  recently  developed  bubble 
breaking  system  described  In  Section  7. 

6.2.1  Four  Atmosphere  Results 

A  total  of  twenty  sets  of  holographic  Interferograms  were  obtained 
of  different  events  at  a  tank  pressure  of  four  atmospheres.  This  tank 
pressure  and  the  Initial  bubble  diameter  of  one  Inch  corresponds  to  a 
Reynolds  number  of  5500.  Two  such  sets  of  Interferograms  are  shown  In 
Figures  12  and  13  for  comparison  purposes.  The  dark  protuberance  in  the 
vertical  views  Is  the  bubble  tube  below  the  vortex.  The  two  vortices  are 
approximately  3.5  and  4  Inches  in  diameter,  and  both  have  risen  about 
four  Initial  bubble  diameters.  It  Is  quite  apparent  that  In  general,  the 
vortices  are  not  Identical  In  shape  or  symmetric  In  a  given  event.  The 
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horizontal  view  of  the  vortex  In  Figure  12  Is  one  of  the  more  symmetric 
In  appearance  while  that  In  Figure  13  Is  one  of  the  more  asymmetric. 

Radial  fringe  number  profiles  for  four  typical  vortex  halves  arc 
shown  In  Figure  14  along  with  a  mean  distribution  based  on  seven  selected 
tests.  The  mean  radial  distribution  of  fringe  numbers  Is  depicted  by 
the  bold  curve  In  the  figure.  The  centers  of  the  vortices  were  chosen 
to  correspond  to  the  minimum  Internal  fringe,  not  to  a  tank  fixed 
coordinate  directly  above  the  bubble  tube.  Presentation  of  the  data  In 
this  manner  tends  to  eliminate  the  effects  of  the  larger  scale  motion 
(l.o..  translational  motion  of  the  vortices). 

A  mean  helium  concentration  profile  for  the  4  atmosphere  data  Is 
presented  In  Figure  15.  This  profile  was  calculated  from  the  mean  fringe 
number  distribution  using  the  Abel  Inversion  program  described  In 
Section  5.1.2.  The  profile  shows  a  helium  concentration  rapidly  Increasing 
from  zero  at  the  vortex  edge  (r/rQ  -  1)  to  a  peak  of  about  1.3  percent  at 
0.75  of  the  normalized  radius  followed  by  a  decrease  to  zero  at  about 
0.25  of  the  vortex  radius.  The  helium  .onccntratlon  becomes  negative 
near  the  center  of  the  toroid  due  to  the  assumption  of  axial  symmetry 
Inherent  In  the  data  reduction  program.*  When  sufficient  test  results  arc 
acquired  and  averaged,  they  will  produce  a  fringe  distribution  corresponding 
to  a  symmetric  vortex  and  consequently  will  be  compatible  with  the  Abel 
Inversion  procedure. 

6.2.2  Eight  Atmosphere  Results 

Ten  sets  of  Interferograms  were  acquired  at  a  tank  pressure  of  eight 
atmospheres.  The  Reynolds  number  for  these  events  was  11,000  based  on 
Initial  helium  bubble  diameters  of  one  Inch.  One  of  the  more  symmetric 
events  Is  shown  In  Figure  16.  This  vortex  Is  approximately  three  Inches 

*)t  must  be  remembered  from  the  governing  equation  for  fringe  shift  that 
fringes  In  Interferograms  correspond  to  changes  In  Index  of  refraction 
Integrated  along  the  path  length  of  the  light  ray  as  It  passes  through  the 
disturbance.  The  assumed  axial  symmetry  Invokes  a  prescribed  path  length 
dependent  only  on  the  diameter  of  the  vortex  and  the  radial  position  at 
which  the  helium  concentration  Is  to  be  calculated. 
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In  diameter  and  has  risen  about  five  Initial  bubble  diameters.  It  Is 
obvious  by  comparing  this  event  with  the  lower  Reynolds  number  Interferograms 
In  Figures  12  and  13  that  the  higher  Reynolds  number  fireball  contains  a 
greater  amount  of  small  scale  turbulent  structure.  In  addition,  the 
fringes  are  closer  together  In  the  higher  pressure  events  which  combines  to 
make  the  data  reading  and  analysis  far  more  difficult  than  In  the  lower 
pressure  cases.  With  the  present  8  atmosphere  data,  It  Is  possible  to 
number  only  the  outermost  fringes.  Preliminary  results  from  a  mean  of 
these  outer  fringe  distributions  were  used  to  calculate  a  peak  helium 
concentration  level  of  slightly  over  one  percent  for  the  8  atmosphere  test 
condition.  Improvements  In  the  holographic  reconstruction  technique  and 
possibly  the  use  of  hydrogen  as  a  replacement  for  helium  In  the  experiments 
are  currently  being  studied  as  means  of  Improving  the  fringe  resolution 
for  the  high  pressure  tests. 

It  Is  Important  to  realize  that  fringe  shifts  In  an  Intcrferogram 
of  a  turbulent  event  result  from  a  combination  of  both  a  mean  and  fluctuating 
density  field  Integrated  along  the  light  pat1".  As  a  consequence,  a  large 
sample  of  data  obtained  at  the  same  test  ccrwiltlons  must  be  ensemble  averaged 
so  that  smaller  scale  asymmetries  due  to  turbulence  will  be  averaged  out. 

A  valid  mean  density  profile  can  then  be  calculated  with  the  Abel 
Inversion  program  since  the  Inherent  requirement  of  axial  symmetry  will 
be  satisfied.  Statistical  studies  are  now  In  progress  to  establish  the 
number  of  data  samples  necessary  to  achieve  desired  levels  of  accuracy 
In  the  reduced  helium  concentration  profiles.  It  Is  now  quite  apparent 
that  50  to  100  or  more  sets  of  repeat  data  will  be  required  at  each  test 
condition  of  Interest. 
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7.  CURRENT  STATUS  AND  FUTURE  PLANS 

A  recent  major  effort  has  been  directed  towards  Improving  the  bubble 
breaking  technique  thought  to  be  the  primary  cause  of  the  larger  asymmetries 
In  the  simulated  fireballs.  Several  different  electrical  discharge  and 
mechanical  devices  wer2  tested  Including  the  use  of  multiple  electrodes, 
but  none  produced  satisfactory  results.  A  mechanical  breaking  technique  of 
dropping  a  small  metal  sphere  through  the  top  of  the  bubble  has  repeatedly 
produced  a  symmetric  helium  release  In  tests  at  one  atmosphere. 

An  Interferogram  obtained  using  this  technique  Is  presented  In 
Flgu  e  17.  This  Interferogram  was  obtained  several  milliseconds  after 
the  bubble  was  broken  ard  shows  the  sphere  just  above  the  bubble  tube. 

As  a  basis  of  comparison,  a  typical  early  time  interferogram  obtained 
using  a  single  electrode  breaker  Is  shown  In  Figure  18.  The  Improvement 
obtained  with  the  new  bubble  bursting  technique  is  obvious.  Two  additional 
Interferograms  that  were  obtained  with  this  technique  were  presented  earlier 
In  Section  4  as  Figures  7  and  8. 

The  prototype  equipment  to  chock  this  bubble  bursting  concept  used 
125  mil  diameter  spheres  and  did  not  allow  the  test,  tank  to  be  pressurized. 
Equipment  Is  currently  being  fabricated  for  the  remote  feed  and  release  of 
smaller  40  mil  diameter  steel  spheres  from  the  top  of  the  test  tank.  Once 
this  new  equipment  Is  Installed  and  tested,  over  50  repeat  sets  of  both 
Interferometry  and  particle  tracking  data  will  be  acquired  at  the  4  atmosphere 
test  condition  described  In  Section  6.2.1.  Repeat  sets  of  similar  data 
will  then  be  obtained  at  8  and  10  atmospheres  pressure  and  finally  at 
10  atmospheres  pressure  In  anSF6  environment  to  provide  an  additional 
increase  In  Reynolds  number. 

Atmospheric  simulation  tests  will  also  be  performed  by  rapidly 
reducing  the  tank  pressure  with  the  vent  system  while  the  torus  Is 
formed  and  rises.  In  this  phase  of  the  experiments  the  effect  of  fireball 
expansion  resulting  from  pressure  variation  with  altitude  will  be  simulated. 
However,  the  effect  of  stratification  or  atmospheric  scale  height  will  not 
be  simulated  in  the  currently  planned  experiments. 
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The  emphasis  in  the  data  reduction  and  analysis  will  initially  be 
oriented  towards  ensemble  averaged  data  and  determination  of  fireball 
mean  velocity  and  density  characteristics.  With  this  data  base 
established,  attention  will  then  be  focused  on  the  primary  subject  of  the 
study,  a  determination  of  the  fireball  turbulent  field  and  its  effect  on 
entrainment  and  mixing. 
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Figure  1.  High  Pressure  Test  Tank  and  Vent  System 
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Figure  3.  Schematic  Diagram  of  Holocamera 
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Figure  5.  Infinite  Fringe  Interferogram  of  Helium  Filled  Soap  Bubble 


Figure  6.  Finite  Fringe  Interferogram  of  Helium  Filled  Soap  Bubble 
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Figure  7.  Infinite  Fringe  Interferogram  of  Bubble 
Shortly  after  Bursting 


Figure  8.  Finite  Fringe  Interferogram  of  Hubble 
Shortly  after  Bursting 
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Figure  9.  Shadowgraph  of  Laminar  Vortex 


Figure  10.  Shadowgraph  of  Turbulent  Vortex 
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Figure  11.  Shadowgraph  of  Turbulent  Vortex  Evolution 
from  Helium  Bubble 
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Horizontal  View 

Figure  13.  Infinite  Fringe  Interferogram  of  4  Atmosphere  Vortex 
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Figure  14.  Radial  Fringe  Number  Distributions  -  Typical  4 
Atmosphere  Events  and  Average 
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Figure  15.  Radial  Helium  Concentration  Profile  -  Averaged  4 
Atmosphere  Data 
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Vertical  View 


Horizortal  View 


Figure  16.  Infinite  Fringe  Interferogram  of  8  Atmosphere  Vortex 
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Figure  17.  Interferogram  of  Helium  Release  -  Figure  18.  Interferogram  of  Helium  Release 
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